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1. Introduction 

The classical integrable systems has been formulated in terms of the classical r- 
matrix In one sense, the system is proved to be integrable when we can perform 
the separation of variables; the reduction of a multi-dimensional system to a set of one- 
dimensional systems (see, for review, Ref. 0). Although the separation of variables 
has been widely known as the Hamiltonian-Jacobi equation, Sklyanin proposed a new 
technique (functional Bethe ansatz), which is closely related with the (quantum) inverse 
scattering method. The functional Bethe ansatz method was first applied to classical 
top 0, and further applications, for the Toda lattice, Gaudin magnet, and Heisenberg 
spin chain, have been done. This technique is a new tool in studying integrable systems. 

We briefly review the separation of variables for the Gaudin magnet in classical 
case. The fundamental formulation of integrable system is based on the classical r- 
matrix structure, 

{L{u),i{v)} = [r{u-v),Liu)+i{v)]. (1.1) 

1 2 

Here we have used the standard notation, L(m) = L(m) ® 1 and L{v) = 1 L('y)- As an 
example of the L-matrix satisfying the linear Poisson structure (|1 . 1|) , we can take 

N c 

LH = z + E— V- (1-2) 

j=i " 

Here (j = 1, . . . , A^) is a classical SL(?t,) spin matrix, and its elements 5*"* (a, b = 
1, . . . ,n) satisfy the Poisson relation, 

{Sf, } = 5jk ■ {S'' Sf - 5'"" Sf) , (1.3) 

and TrS| is the Casimir element. We suppose that matrix Z is traceless, TrZ = 0. 
The L-matrix ( |1.2| ) appears as a quasi-classical limit of the inhomogeneous Heisenberg 
XXX-spin chain [§,0, and satisfies the linear Poisson relation ( |1 . 1| ) with the classical 
r-matrix, 

P 

r{u) = — . (1.4) 
u 

Matrix P means a permutation matrix, p°'^''^'^ = S°''^S'"^, which satisfies Px^y = y®x. 
Note that r-matrix ( [L.4|) is a rational solution of the classical Yang-Baxter equation 0, 

[r23(t^), ri2(M)] + [r23(w), ri3{u + v)] + [ris{u + v), ruiu)] = 0. (1.5) 

Once the L-matrix satisfies the Poisson structure ( p. . 1|) , the model can be proved to 
be integrable in the Liouville's sense. If functions of the L-matrix, Tm{u), are defined as 

TUu) = -TTL{ur, (1.6) 

m 
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one can see that Tm{u)^s are spectral invariant, 

{ti{u),tM} = 0. (1.7) 

The identity ( p..7| ) shows that rm(M)'s are generating functions of the constants of motion. 
The first non-trivial invariant follows from T2{u); 

Mu)^l^Z'.E^^.l^^^^, (1.8) 

where Hj is the hamiltonian of the SL(?2) Gaudin magnet, 



Hj = Tr ZSj + J2 —■ (1-9) 

k^j ~ 

This model was introduced by Gaudin as an integrable spin system with long-range 
interaction p[. Due to the involutiveness of Tm{u), one sees that the hamiltonian of the 
SL(r2) Gaudin magnet is Poisson commutative, 

{H^,Hk}=0, forj,A; = l,...,iV. (1.10) 

The complete integrability of the model in Liouville's sense can be proved directly 
from (|1.7|); when we introduce quantities t^j by 

1 N m _a 

Tm{u) = - Tr Z- + 5: 5: (1.11) 

j=l a=l \ J y 

we can see that the quantities t^j {m = 2, . . . , N; j = 1, . . . ,n; a = 1, . . . ,m — 1) form 
a commutative family of Nn{n — l)/2 independent hamiltonians. 

For this type of the Gaudin magnet (|1.9| ), the separation of variable (functional 
Bethe ansatz) has been widely studied in both classical and quantum cases P,p|,p^,pA|, 
r2|,|13],0. Let ^(L) and B(L) be certain polynomials of degree n{n — l)/2 in matrix 



elements Lab- When we define variables Xj and pj by 

B{L{x,)) = 0, p,=A{-L{xj)), (1.12) 



one sees that variables Xj and pj are canonically Poisson conjugate [pr5|,p!6|,p!7 



{xj,Xk} = 0, {Pj,Pk} = 0, {Pj,Xk} = 6jk. (1.13) 

This analysis, which is called the separation of variables, makes it possible to calculate 
the energy spectrum for the quantum Gaudin magnet. 

By this way, we can perform the separation of variables for integrable systems 



formulated in the linear Poisson relation (1.1). Some of integrable systems, e.g 



nonlinear integrable equations on finite segment, are formulated in terms of another 
Poisson structure; there exists 'rs'-Poisson structure |]18|,0, 



{t{u),i{v)} = [r{u-v),t{u)+L{v)] + [s{u + v), -i{u)+ i{v)].{l.U) 
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Here r and s are matrix structure constants. Remark that s-matrix depends on sum of 
spectral parameters while r-matrix on difference. The Poisson structure ( |1.14|) can be 
viewed as a classical limit of the boundary Yang-Baxter equation , which is used to 
formulate the quantum spin chain with open boundary. 

In this paper we shall study the separation of variables for BC-type integrable system 
formulated by the rs-Poisson structure ( |1.14| ). In section 2, we introduce the BC-type 
SL(n) Gaudin magnet. We give classical r- and s-matrices, and prove the integrability. 
We relate the hamiltonian of the BC-type SL(2) Gaudin magnet with the generalized 
Neumann model in section 3. The separation of variables is also studied. In section 4, we 
pay attention on quantum case. The energy spectrum is given based on the separation 
of variables. Section 5 is devoted to conclusion and discussion. 



2. Gaudin Magnet with Boundary 

The quantum Gaudin magnet, whose hamiltonian has a form ( |1.9| ), was first introduced 
in Ref. as an integrable spin system with long-range interaction, and solved by 
use of the coordinate Bethe ansatz. As was reviewed in the previous section, this 
original Gaudin magnet can be formulated with the linear Poisson structure in classical 
case ( p..!] ). In this section, we consider the SL(n) Gaudin magnet with boundary (BC- 



type Gaudin magnet). The dynamical variables of this model are the classical SL(?2) 
spin Sj^ {a,b = 1, . . .n; j = 1, . . . ,N) satisfying the Poisson bracket ([T7^ ). Consider the 



modified L-matrix, 

uu)=j:(-^+-^), (2.1) 

where the "reflected" classical spin S is defined as 

■^ab ^ ^_y+bgab_ ^2.2) 



See the difference from the usual L-operator ( |1.2| ). The second term in (2J.) is due to 
the effect of the reflection; the classical spin Sj is located at coordinate Zj, while the 
"reflected" spin Sj is at —Zj. For this reason we say that the system has "boundary". 
One can check easily that the modified L-operator ( |2.1| ) satisfies the linear Poisson 
structure ( |1.14| ) with the rational r-matrix ( |1.4| ) and s-matrix. 



s{u) = -, (2.3) 
u 

where we use the notation, p"'''''^'^ = ^_Y+b§o-d^bc^ 

Let us define functions Tm{u) of the matrix L(m) ( |2.1|) as. 



TUu) = -TTL{ur. (2.4) 
m 
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From the Poisson structure ( |1.14| ), it is shown that functions Tm{u) are the spectral 
invariants of dynamical system, 

{ti{u),t^{v)} = Q. (2.5) 

The first nontrivial invariant is given from T2{u) as, 

- £ ... '/w!'^ + \ £(7^ + (2.6) 



\u - Zj){u + Zj) 2 V(m - Zjf {u + zj)^ 



where Hj has a form. 



Tr(S,S,) ^ £./ Tr(S,S.) ^ Tr(S,S.) Y ^^.7) 

2 V Zj Z]^ Zj ~\- Zi^ J 

The involutiveness of the spectral invariants Tm{u) indicates the Poisson commutativity 
of the hamiltonian Hj , 

{H,,Hk] = Q, forj,A; = l,2,...,iV. (2.8) 

We call Hj the hamiltonian of the SL(n) Gaudin magnet with boundary. Different from 
the original Gaudin magnet (|1.9| ), the hamiltonian Hj includes interacting terms between 
classical spins Sj and "reflected" spins Sj. By calculating the residues of Tm{u) ( |2.4|) as 
in the case of original Gaudin magnet (|1 . 1 1 ) , one can get "higher-order" hamiltonian 
r°j- of the Gaudin magnet. The existence of commutative family {t^j} supports the 
complete integrability of model in Liouville's sense. 



3. Separation of Variables 

In this section we first show the isomorphism of the iV-site SL(2) Gaudin magnet and the 
iV-dimensional generalized Neumann model, and then study the separation of variables 
based on the technique of Sklyanin. In SL(2) case we can define the classical spin 
matrices Sj and Sj in the L-operator ( p.l| ) as, 

S,= (|( fj. ), S, = a^S,a^ (3.1) 

where is denoted as Pauli spin matrix. These spin variables satisfy the following 
Poisson relations; 

{Sj, S^} = ±5jk Sf, {57, "S*^ } = 25jfc Sj. (3.2) 

Above Poisson structures for spin variables can be realized with new variables Xj and 

Pj as. 
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where {xj,pj\j = 1,...,A^} are canonical variables satisfying the Poisson relations; 
{xj, Xk} = {Pj,Pk} = 0, and {xj,pk} = Sjk- In this case the Casimir element is set to be 
zero, TrSj = 0. In terms of canonical variables, one can obtain the hamiltonian from 
the spectral invariant T2{u), 

^ z- 

= E7 V7 ^H,, (3.4) 

where Hj is calculated as, 

jj _ 1 ^( {Pj^k - XjPkf jpjXk + XjPkf \ ^^^^ 

Zj 2 \ Zj Zf^ Zj -\- Z]^ J 

This hamiltonian can be viewed as generalization of the Neumann model |]^,|n|. This 
proves the fact that the A^-dimensional generalized Neumann model ( |3.5D is isomorphic 
to the A^-site BC-type SL(2) Gaudin magnet (^). 

Now we study the separation of variables for the generalized Neumann model ( |3.5D . 
Define A{u) and B{u) as 

A(n)=LnH, B{u) = L^2{u). (3.6) 



The linear Poisson structure ( |1.14|) includes relations between functions A{u) and B{u) 



{A{u),A{v)} = Q, {B{u),B{v)] = Q, 

2v / \ (3-7) 

{A{u),Biv)} = ^^—-{b{u) - B{v)). 

[U — V)[U + V) ^ ' 

Both matrix elements, A[u) and B{u), Poisson-commute among themselves. 

We choose separable coordinates as zeros of the off-diagonal element of L-operator, 



B{±Ua) = 0, for a = l,2,...,iV- 1. (3.8) 

Note that u = —Ua is solution ii u = Ua solves equation B{u) = 0. By use of a set of 
variables Ua, we further introduce the canonical variables by 

Va = A{u^). (3.9) 

From the Poisson relations ( p.7|) one can see that variables and are canonically 
conjugate, 

{Ua,Up} = 0, {Va,Vf3} = 0, {Ua,Vp}=Sap. (3.10) 

The first two Poisson relations are easy to be proved. The third relation follows from 
{uo„vp}= lim } ( — ^ \ — ]b{u). 



u^ua B' {u)\Up — U Up + u 



The Poisson relations ( p.lO| ) show that Ua and Va are canonically conjugate variables, and 
that the generalized Neumann model (|3.5| ) is separated by transforming the dynamical 
variables as 

{^j^Pj\j = l,---,^} ^ {Ua,Va\j = l,...,iV- 1}. 

With variables Ua and Va, the action W of the generalized Neumann model (|3.5| ) is 
written as separated form, 

N-l „ 

W=Y. Vo^duo.. (3.11) 

The separated variables Ua and Va can be written explicitly in terms of Xj and pj. By 
definition ( p.8D we can solve the coordinates pj as 



N-l 



P? = ^ ■ ^ , (3.12) 



with z = J2j Zjp'j- The canonically conjugate variables Va are also solved as 

N 

^ AM = - E 7 (3.13) 

j^l [Ua - Zj){Ua + Zj) 

Remark that the relation between {ua,Va} and the hamiltonian Hj is given from the 
spectral invariants T2{u = Ua), 



V. 



f^l - Zj){Ua + Zj) 

4. Quantum Case 

In this section we consider to quantize the BC-type Gaudin magnet ( p.7|) . For brevity 
we study SL(2) case. We set the L-operator as 

uu) = j:(-^+-^), (4.1) 



J-- 



"[XU — Zj U + Zj 



where the spin operator Sj and reflected spin operator Sj are defined as 

«^=(| (4.2) 



Here operators S'| and 5*^^ denote bases of the su(2) Lie algebra, and they satisfy the 
commutation relations, 

[Sj, S^] = ±Sf6jk, [Sj', Sf. ] = 2Sj6jk, 

1 (4.3) 
{S;Y + -iS^S- + S-SD = + l), e Z+/2. 

We have set ij as spin of the j-th site. One can check from direct calculations that the 
L-operator ( [4. 1]) satisfies the quantum analogue of the linear Poisson structure (|1.14| ), 

[iiu),tiv)] = [r{u - v),i{u)+ tiv)] + [s{u + v), -i{u)+ i{v)]. (4.4) 

In this case the constant r- and s-matrices are defined as 

P 1 1 

r(u) = , s(u) r(u) a"" . (4.5) 

u 

The conserved operators are generated by the same way in the classical case; trace 
of the L-operator (|2.4|) . We get the first non-trivial operator from 7^2 (u), 

f2(«) = ^TrL(w)2 

Here the quantum operator Hj is the hamiltonian of the quantum BC-type Gaudin 
magnet, 

4. = ^ + f: C^^^ + ^^^^) ■ (4.7) 

2 Zj ^ / ■ \ Zj Zk Zj ~\- Zk J 



From the commutativity of the generating function f2{u), we can see that the operators 
Hj are commutative, 

[H„Hk]=0, ioTj,k = l,...,N, (4.8) 

which proves the quantum integrability of the BC-type Gaudin magnet. Note that 
the operator Hj ([4.7|) has been appeared in recent studies of the generalized Knizhnik- 
Zamolodchikov (KZ) equation pli|2^ . 

The separation of variables for quantum case can be performed as follows. When 
we define operator A{u) and B{u) as 

A{u)=Lniu), B{u) = Li2{u), 

we obtain from the quantum rs-structure ( ^.4|) that the commutation relations among 
operators A{u) and B{u) can be written as, 

\A(u),A(v)] = 0, |B(«.),B(k)] =0, 
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All the calculation is essentially same with the classical case. In the quantum case, we 
can also introduce the "canonical operators" Ua and Va by, 

Bi±Ua) = 0, 

(4.10) 

These operators and Va satisfy the commutation relations, 

[Ua,Up] = 0, [Va,Vi3] = 0, [Ua,Vp] = 5af3. (4.11) 

To perform the separation of variables ( [4.10 ) for the quantum BC-type Gaudin 
magnet, we use the realization of the spin operators, and Sf, as 

With this realization, the functional equation, B{±Ua) = 0, does not include differential 
operator, and can be solved easily. The result is 

Af-l 

n -Ua){zj +ua) 

= ^^^^ , (4.13) 

where we set z = J2kZkXk- When we change the variables from {xj} to we can 

see that operator Va is represented in terms of 

Va = -l^+MUa), (4.14) 
OUa 

where we use function A{u) as 

A{u) = y2- ^ -. (4.15) 

Remark that identification of operator Va ( |4.14] ) is consistent with the commutation 
relations ( [4.11| ). A set of operators {ua,Va} is called the separated operator. 

We have completed in separating variables for quantum SL(2) Gaudin magnet. In 
the rest of this section, we show that the energy spectrum for the Gaudin magnet can be 
calculated from the functional Bethe ansatz. By definition of the generating function, 
^2(w) = |TrL(M)2, one can see that 

vl-f2M = 0, (4.16) 

which corresponds, in classical case, to ( p.l4 ). With the operator realization of Va 
obtained in ( |4.14| ), we can read off the identity ( [4.16| ) as the differential operator for the 
separated spectral problem, 

ij"{u) - 2A(m)V^'(m) + (a2(m) - A'{u))-^{u) = T2{u)iIj{u), (4.17) 
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where T2{u) is an eigenvalue of operator f2{u) (4.6). This equation can be seen as a 
generahzed Lame equation [10|. To solve this second order differential equation ( |4.17|) , 
we assume that the wave function ipi'^) is a polynomial of u, and that zeros of ip{u) are 
denoted as ±Xa M, 



M 



(4.18) 



a=l 



Substituting the wave function ■?/'(«) into differential equation ([4.17|) , we can see that 
the eigenvalue Ej of Hj is given by, 



N 



4z 



(4.19) 



where function x[u) is defined from the wave function ^/'(m) as 



d-u 



log '?/'(«) 



M 

E 



2u 



(4.20) 



Notice that the zeros of the wave function tp{u) should be fixed to satisfy a set of 
equations, 



A(A, 



M 

E 



2 Aq, 



+ 



1 



(A«-A^)(A« + A^) 2 A 



for a = 1, 



(4.21) 



which follows from the conditions in cancellation of the residues at m = Aq, in (|4.17| ). This 
equation is a quasi-classical limit of the Bethe ansatz equation for open boundary spin 
chain, and plays a crucial role in construction of the integral solution for the generalized 
KZ equation p2| . 



5. Discussion 



We have introduced the generalized Gaudin magnet in this paper. The hamiltonian is 
written as 

Tr(S,S,) , ^^ Tr(S,S,) , Tr(SA) ^ 

This model can be regarded as the Gaudin magnet with boundary, and be formulated 
in terms of "classical" reflection equation (|1.14|) . As in the case of the original Gaudin 
magnet, this model has many interesting aspects; in particular for SL(2) case, the model 
is proved to be isomorphic to the generalized Neumann model. In both classical and 
quantum cases we have performed the separation of variables. In this analysis the 
eigenvalue problem can be reduced to the second-order differential equation (Lame 
equation). We can obtain the so-called quasi-classical Bethe ansatz equation from this 
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differential equation. The XXZ-Gaudin magnet with boundary will be analysed in the 
same method. 

The point of the separation of variables (functional Bethe ansatz) is to take zeros 
of the wave-function ip{u)] the zeros may be identified with the "rapidities" of the spin- 
wave from the view point of the inverse scattering method. This kind of analysis was 



used in recent studies of the Asbel-Hofstadter problem [Q. From the viewpoint of 
the g-polynomial theory the Askey- Wilson polynomial may be related with quantum 
rs-structure fM . 
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